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ABSTRACT 

In a medium where a photoreactive species produces N equally pho- 
toreactive radicals, by combining the Lambert-Beer law for the decay 
of light intensity, as a beam propagates into the medium, with the 
kinetic equation, one obtains a generalization of the Wegscheider 
equation from photobleaching theory. It is shown that this equation 
can be solved exactly, but implicitly, and can be reduced to a first 
order ordinary differential equation in a single reduced variable. 
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Introduction 



As a light beam propagates into a medium containing a various photoreactive 
species the beam intensity decreases with depth due to absorption proportional 
to the concentration of the active molecules as is described by the Lambcrt-Bccr 
law [1] . By combining this law with the kinetic equation one obtains a non- linear 
first order integro-partial differential equation. Situations where the N radicals 
released lead to desirable effects, such as bleaching and polymerization, have 
been studied for many years. An extensive list of references is presented in [2] . 
For the case of a single species the Lambert-Beer law for the light intensity 
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I{x,t) = Ioexp[—a / C{u,t)du] (1) 
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baI{x,t)C{x,t) (2) 



combined with the kinetic equation 

dC{x,t) _ 
dt ~ 

gives Wegscheider's equation [3] 

dC(x,t) 

— — — = — ^a/oexp[— a / C{u,t)du]C{x,t) 
(^t Jo 

C{x,0) = Co. (3) 

whose solution has been known for many years [2] . Here C is the concentration, 

a is the absorption coefficient, and (j) the quantum yield. For the case of a 
slab < a; < L, in terms of reduced variables T = (j)at, z = xjL, Sq{z,T) = 
C{x, t)/Co, 7 = aCoL, the solution is 

So(z.T)^[l + e-'^-ie^ -1)]-^ (4) 

which is examined thoroughly in [2]. 

Systems are presently under investigation [4] where, not only the initial 
species, whose concentration is C, but its N reaction products all have the same 
absorption coefficient. A similar analysis leads to the more general equation, 
expressed in terms of the dimensionless variables given above. 



dC{z,T) 
df 



exp[7 [ C{u,T)du]C{z,T) 
Jo 



C(z,0) = l. (5) 
The aim of this note is to present the exact, albeit implicit, solution to (5). 
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Calculation 

We first introduce the cumulative concentration a{z,T) = S{u,T)du for 
which a{z, 0) = ^ and a{0, T) = 0. Then by integrating both sides of (5) over z 
and then differentiating with respect to z, one finds 

[a{z, T) - Nz] = _e-T['^(^'^)-^"l |- [(t(z, T) - Nz] . (6) 



dzdT' ^ ' ' ' dz 
That is, for f{z,T) = a{z,T) — Nz one has the partial differential equation 

^2 



X 



rp 1- Ns-e- 



or, since S = e ^ , following a simple change of integration variable, 

-T 

du 



lsiz,T) + A^log u-u) 
Finally, in terms of the new variables 



du 



(7) 



with f{z, 0) = -{N - l)z, /(O, T) = 0. 

Next, we introduce V{x, T) = — log[/z + N] to obtain 

which, after integration with respect to T becomes 

^^=,(l-W-e-) 
oz 

V{0,T)=T. (9) 
From (9) we get the implicit relation 

^^(^'^) ds 



- 70 = (10) 



7^ = 0. (11) 



/In 2 1 - - 

^ = JZ-T (12) 

we find that the solution to (5) has the implicit representation 

du ^ ^ 

J 1/2 u{l+N\n.u — u) 



(13) 
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Discussion 



Let us first look at the case = 0, which renders (5) equivalent (under 
S = —So) (mathematically, but not physically) to the Wegscheider equation [1]. 

By explicit integration t = ln(e-^ — 1), S'(^) = e^(e^ + 1)"-'^ and wc recover the 
solution (4). Note that (13) is equivalent to the first order ordinary differential 
equation 

5'(0=C(l + A^ln^-0 (14) 

subject to an appropriate initial condition. Eq. (14) should be useful for obtain- 
ing series approximations to S{z, T). It also indicates that for > S{z, T) is 
nonanalj^ic along the trajectory = r. 
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